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ABSTRACT 

Lensing probabilities of quasars with image separations greater than A6 
and flux density ratios less than q T are calculated by foreground dark mat- 
ter halos in a flat, cosmological constant dominated (ACDM) universe . The 
mass density of the lenses is taken to be the Navarro- Frenk- White (NFW) 
profile on all mass scales, plus a central point mass for low-mass halos with 
M < M c = 5 x 10 13 h~ l M Q . We introduce a quantity M e g, which is a point mass 
ranging from 1 to 1000 times the mass M. of a supermassive black hole (SMBH) 
inhabiting the center of each galaxy, to describe the contributions of galactic 
central SMBHs and galactic bulges to lensing probabilities. The lensing cross 
section and thus the lensing probability are quite sensitive to the flux density 
ratio q T of multiple images in our calculations. It is shown that, to reproduce 
the lensing survey results of JVAS/CLASS for q T < 10, about 20% of the bulge 
mass is needed as a point mass for each galaxy. Since there is still considerable 
uncertainty regarding the value of the spectrum normalization parameter erg, we 
investigate the effect of varying this parameter within its entire observational 
range (from 0.7 to 1.1), and find that low o~$ values (< 0.7) are ruled out, and 
the best fit value is as — 1.0. 

Subject headings: cosmology: theory - cosmology: observations - galaxies: bulges 
- galaxy: center - gravitational lensing 



1. Introduction 

Gravitational lensing provides us a powerful probe of the mass distributions of the 
universe. By comparing the lensing probabilities predicted by various cosmological models 
and the density profiles of lenses with observations, we are able to test the mass distribution 
of dark matter (CDM) halos and in particular, the inner density slope in the sense that the 
JVAS/CLASS radio survey (Browne & Meyers 2000; Helbig 2000; Browne et al. 2002; Myers 
et al. 2002) has provided us the observed lensing probabilities at small image separations 
(0.3" < A6 < 3"). 
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It is well known that, Cold Dark Matter (CDM) model has become the standard theory 
of cosmological structure formation. The ACDM variant of CDM with Vt m — 1 — Q,\ !=» 0.3 
appears to be in good agreement with the available data on large scales (Primack 2002). On 
smaller (sub-galactic) scales, however, there seems to be various discrepancies. Issues that 
have arisen on smaller scales have prompted people to propose a wide variety of alterna- 
tives to the standard CDM model, such as warm dark matter (WDM) and self-interacting 
dark matter (SIDM). Now that problems arise from galaxy-size halos and the centers of 
all dark matter halos, high-resolution simulations and observations are the final criterion. 
Recent highest-resolution simulations appear to be consistent with NFW (Klypin, Zhao, & 
Somerville 2002; Power, Navarro, & Jenkins 2002) down to scales smaller than about 1 kpc. 
Meanwhile, a large set of high-resolution optical rotation curves has recently been analyzed 
for low surface brightness (LSB) galaxies, which suggests that the NFW profile is a good fit 
down to about 1 kpc. Although further simulations and observations, including the obser- 
vations of CO rotation curves (Bolatto et al. 2002), may help to clarify the issue, it is likely 
that both WDM and SIDM are probably ruled out, while the predictions of ACDM at small 
scales may be in better agreement with the latest data. 

Unfortunately, if dark halos are modelled with NFW profile on all mass scales, they 
will produce too few small angular separation images in ACDM model as compared with the 
results of JVAS/CLASS (Li & Ostriker 2002). A possible solution to the problem is to modify 
the inner structure of dark halos by introducing baryonic cooling and compression (Porciani 
& Madau 2000; Kochanek & White 2001; Keeton 2001; Sarbu, Rusin, & Ma 2001; Li & 
Ostriker 2002; Oguri 2002). Namely, there exist two types of dark halos: small mass halos 
(galactic size) with a steep inner density slope (singular isothermal sphere, SIS) and very 
massive halos (e.g. galaxy clusters) with a shallow inner slope (NFW). This can reproduce 
the observational data of JVAS/CLASS (Li & Ostriker 2002; Sarbu, Rusin, & Ma 2001). 
One may also consider a unified model in which dark halos are composed of an NFW-like 
component and a bulge component for galaxies as lenses(Chen 2003). If the NFW profile is 
believed to be universal, this model allow us to constrain the structure of the galactic centers 
using the strong lensing surveys. 

In this letter, we investigate the plausibility of the NFW+Bulge model by fitting the 
observational data of JVAS/CLASS in a much more accurate way to improve our previous 
work (Chen 2003). Furthermore, we emphasize the importance of the flux density ratio q r 
of the two images produced by a central point mass in each galaxy for the predicted results. 

In our model, there are two important issues when galactic bulges are involved: First, 
the presence of supermassive black holes (SMBHs) at the centers of most galaxies appears 
by now firmly established (Melia & Falcke 2001, and references therein). SMBH masses are 
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estimated to be in the range 10 6 — 1O 9 M , and are correlated with the masses and luminosities 
of the host spheroids and, more tightly, with the stellar velocity dispersions (Magorrian et al. 
1998; Ferrarese & Merritt 2000; Ravindranath, Ho, & Filippenko 2001; Merritt & Ferrarese 
2001a,b; Wandel 2002; Sarzi et al. 2002). Recent high-resolution observational data give 
M./Mb u ig e ~ 10 -3 (Merritt & Ferrarese 2001c). Furthermore, Ferrarese (2002) finds a relation 
between masses M. of SMBHs and the total gravitational masses of the dark matter halos 
as (M./1O 8 M ) ~ O.O46(M DM /1O 12 M ) L57 . Since these correlations extend well beyond the 
direct dynamical influence of the SMBH it seems likely that there is a close link between the 
formation of both the SMBH and its galaxy (Silk & Rees 1998; Adams, Graff, & Richstone 
2001; Haehnelt & Kauffmann 2001; Islam, Taylor, & Silk 2002; Madau & Rees 2001; Menou, 
Haiman, & Narayanan 2001; Schneider et al. 2002). So we can simply add galactic SMBHs 
as point masses to the NFW density profile when we calculate lensing probability. 

Second, galactic bulges, in which multiple black holes may form and inhabit (e.g., 
Haehnelt & Kauffmann 2002), can also contribute to the lensing probabilities at small image 
separations. Since light rays are affected by the mass within the sphere of their impact 
distances, we can attribute the light deflections induced by such a mass to an effective point 
mass, which is referred to as M e g in this paper. We thus avoid the complexity of gas com- 
pression and mass distributions of galactic bulges. Of course, we are unable to reveal the 
detail structures of the inner part of galaxies at the same time. We argue that this may be an 
adequate way if we want to compare our predictions with the results of strong gravitational 
lensing survey. So, the total central pointlike masses by including central SMBHs and bulges 
range from 1 — 1000M.. The upper limit of M eS is 1000M. because the bulge mass correlates 
linearly with SMBH mass as M./M bu i gc m 10~ 3 . 



2. Lensing equation and probabilities 

We choose the most generally accepted values of the parameters for the ACDM cosmol- 
ogy, in which, with usual symbols, the matter density parameter, vacuum energy density 
parameter and Hubble constant are respectively: tt m = 0.3, = 0.7, h = 0.75. The NFW 
density profile is pnfw = Psft/Vi 7 " + r s) 2 ]> where p s and r s are constants. We can define the 
mass of a halo to be the mass within the virial radius of the halo r v ; r : Mom = 4"7rp s r 3 /(ci), 
where /(ci) = ln(l + Cl ) - Ci/(1 + c 1 ), and Cl = r vir /r s = 9(1 + z)-\M/l.h x lO^/i^Me)" 1 
is the concentration parameter, for which we have used the fitting formula given by Bullock 
et al. (2001). 



The surface mass density for a halo as lens is 

S(f) = M cff 5 2 (f) + S NFW (f), 



(1) 
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where x = \x\ and x = |*/r s , is the position vector in the lens plane. S 2 (x) is the two 
dimensional Dirac-delta function, and £ NFW (:r) is the surface mass density for an NFW 
profile. M c g is a point mass ranging from 1 to 1000 times the mass M. of a supermassive 
black hole (SMBH) inhabiting the center of each galaxy, to describe the contributions of 
galactic central SMBHs and galactic bulges to lensing probabilities. For galaxy cluster, 
M c g = 0. The lensing equation for this model is then 

y = x-fi s , (2) 

x 

where y — \y\, fj — yD^/D^ is the position vector in the source plane, in which D§ and 
are angular-diameter distances from the observer to the source and to the lens, respectively. 
It should be pointed out that, since the surface mass density is circularly symmetric, we 
can extend both x and y to their opposite values in Eq.(2) for convenience. The parameter 
/4 — 4ps r s/S cr is independent of x, in which E cr is critical surface mass density. The term 
/ eff = 2.78 x 10- 4 /(ci)M?5 57 (Meff/M.), where M 15 is the reduced mass of an NFW halo 
defined as M i5 = M D M/(lO 15 h _1 M ), stands for the contribution of a point mass M eff , and, 
of course, / e ff = for cluster-size lenses. The function g(x) stands for the contribution of 
the NFW halo, and it has an analytical expression originally given by Bartelmann (1996). 

When the quasars at redshift z s — 1.5 are lensed by foreground CDM halos of galaxies 
and clusters of galaxies, the lensing probability with image separations larger than A9 and 
flux density ratio less than q r is (Schneider, Ehlers, & Falco 1992) 

f Zs dDi (z) r°° 

P(>A6,<q r )= ft '- dz / n{M,z)a(M,z)B(M,z)dM, (3) 

Jo dz J 

where -Dl(z) is the proper distance from the observer to the lens located at redshift z. The 
physical number density n(M, z) of virialized dark halos of masses between M and M + dM 
is related to the comoving number density n(M, z) by n(M, z) = n(M, z)(l + z) 3 , the latter 
is originally given by Press & Schechter (1974), and the improved version is n(M,z)dM = 
(po/M)f(M,z)dM, where po is the current mean mass density of the universe, and / = 
(0.301/M)(dlnA z /dmM)exp(-| ln(A z /1.68) + 0.64| 3 - 88 ) is the mass function for which we 
use the expression given by Jenkins et al. (2001). In this expression, A z = 5 c (z)/A(M), 
in which 5 c (z) is the overdensity threshold for spherical collapse by redshift z, and A(M) 
is the rms of present variance of the fluctuations in a sphere containing a mean mass M. 
The overdensity threshold is given by 5 c (z) = 1.68/ D(z) for the ACDM cosmology (Navarro, 
Frenk, & White 1997), where D(z) = g[fl(z)]/[g(Q m )(l + z)\ is the linear growth function 
of the density perturbation (Carroll, Press, & Turner 1992), in which g(x) = 0.5x(l/70 + 
209x/140-a; 2 /140+a; 4 / 7 )- 1 and Q(z) = tt m (l + z) 3 /{l-tt m + Q m (l + z) 3 ]. When we calculate 
the variance of the fluctuations A 2 (M), we use the fitting formulae for CDM power spectrum 
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P(k) = AkT 2 (k) given by Eisenstein & Hu (1999), where A is the amplitude normalized to 
o" 8 = A(r M = 8/i -1 Mpc) given by observations. Note that the mass of an NFW halo is taken 
to be the virial mass M DM = 47r<5 vil .pr^ ir /3, where 5 vir p is the average density within r vir . 

The key step in working out the final results of lensing probabilities is how to calculate 
the lensing cross section cr(M, z) in Eq. (3). Since we are interested in the lensing probabili- 
ties with image separations larger than a certain value AO (ranging from ~ 10 arcseconds, 
for example) and flux density ratio less than q T , the cross section should be defined under two 
conditions. The first condition is used to define the cross section of cluster-size NFW lenses, 
for which, multiple images can be produced only if \y\ < y cr , where y cr is the maximum value 
of y when x < 0, which is determined by dy/dx = when f c s = in Eq.(2). On the other 
hand, because of the existence of the central point mass, theoretically, every galaxy-size lens 
will always produce two images. So the first condition fails in this case and we need the sec- 
ond condition to define the cross section of galaxy-size lenses, which is the allowed upper limit 
of flux density ratio of lensing images in any lensing survey experiments. So, the flux density 
ratio q r for the two images is just the ratio of the corresponding absolute values of magnifica- 
tions (Schneider, Ehlers, & Falco 1992; Wu 1996), q r = where n + [y(x)} = (f^) a . >0 
and H-[y{x)\ = (| J) x<0 - So y CT is determined by |/x+(y cr )| = <? r |/x_(y cr )|. The cross sec- 
tion for images with a separation greater than AO and a flux density ratio less than q r is 
(Schneider, Ehlers, & Falco 1992) 

( y 2 cr , for AO < A0 ; 

a(M, z) = nr 2 J(M - M min ) x 1 y 2 CI - y 2 A9 , for A0 < AO < A0 yci ; (4) 

U for A0 > A0, cr , 

where is a step function, and M min is the minimum mass of halos above which lenses 
can produce images with separations greater than AO. From Eq.(2), an image separation for 
any y can be expressed as A6{y) = r s Ax(y)/D^, where Ax(y) is the image separation in lens 
plane for a given y. So in Eq.(4), the source position y^g, at which a lens produces the image 
separation AO, is the reverse of this expression. And A#o = A0(0) is the separation of the two 
images which are just on the Einstein ring; A0 yci = A0(y CT ) is the upper-limit of the separa- 
tion above which the flux ratio of the two images will be greater than q T . Note that since M DM 
(Mis) is related to A# through r s = (1.626/c 1 )(M 1 1 5 /3 / [Q m (l + zf + fi A ] 1/3 )/ i - 1 Mpc (Li & 
Ostriker 2002), we can formally write M D m = M D M(A0(y)), and determine M min for galaxy- 
size lenses by M min = M D M(A^(y cr )), and for cluster-size lenses by M min = M D m(A^(0)). 
In the latter case we have used the fact that the separation of the outermost images is 
insensitive to the value of y in cluster-size NFW lenses. 

To compare the predicted lensing probabilities with the combined data from JVAS/CLASS, 
magnification bias must be considered. For the JVAS/CLASS sample, we use the result given 
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by Li & Ostriker (2002): B(M,z) = 2.22A 1 m 1 (M, z), where A m (M, z) = Ax(y = 0)/y cr . 

One major uncertainty in the estimate of P(> AO, < q T ) by the NFW halo arises from 
the assignment of the concentration parameter C\ to each halo of mass M. There exist 
several empirical fitting formulae or analytic models to fulfill the task. However, for a given 
halo mass and redshift, there is a scatter in c\ = r vir /r s value that is consistent with a log- 
normal distribution with standard deviation a c = A(logc) ~ 0.18 (Jing 2000; Bullock et al. 
2001). We take into account the scatter in c\ by averaging the lensing probability with the 
log-normal distribution (the right panel of figure 1 and all the panels of figure 2) 

Another Major uncertainty in predicting P(> AO, < q T ) arises from the considerable un- 
certainty regarding the value of the CDM power spectrum normalization parameter er 8 , so it 
would be useful to consider the effect of varying this parameter within its entire observational 
range, roughly <r 8 = 0.7 ~ 1.1 (see Fig. 2). 

3. Discussion and conclusions 

The lensing probabilities predicted by Eq. (3) and calculated from the combined 
JVAS/CLASS data are compared in Fig.l. Since we are interested in the degeneracy between 
q r and M e g in matching the predicted results to observational data, we have calculated the 
lensing probabilities for four different values of q r and their corresponding values of M e s , as 
indicated in Fig. 1. We have assumed a "cooling mass" of M c ps 5 x 10 13 /i _1 M , above which 
the lenses are assigned the NFW profile and below which the lenses are the NFW+point 
mass. The predicted lensing probabilities are quite sensitive to M c , the value of M c used 
here is higher than those preferred by other authors (Kochanek & White 2001; Li & Ostriker 
2002; Sarbu, Rusin, & Ma 2001). A lower value (e.g., M c « 1O 13 /i _1 M ) requires a larger 
M e fj for each q r , however, this will not affect our conclusions, since our main goal is not to 
"measure" the exact value of M e $. The combined JVAS/CLASS survey is now completed, a 
subset of 8,958 sources form a well-defined statistical sample containing 13 multiply-imaged 
sources (lens systems) suitable for analysis of the lens statistics. One of the four obser- 
vational selection criteria for this "well-defined" sample is: the image components in lens 
systems must have separations > 0.3 arcsec and the ratio of the flux densities of the brighter 
to the fainter component in double-image systems must be q T < 10 (Chae et al. 2002). The 
observed lensing probabilities can be easily calculated: P bs(> AO) = N(> A#)/8958, where 
N(> AO) is the number of lenses with separation greater than A# in 13 lenses. P Q bs(> AO) 
is plotted as a histogram in both panels of Fig. 1. 

First of all, as shown in Fig. 1, when averaged over concentration parameter c\ with the 
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Fig. 1. — Predicted lensing probability with image separations > A6 and flux density ratios 
< q r in ACDM cosmology. The cluster-size lens halos are modelled by the NFW profile, 
and galaxy-size lens halos by NFW+BULGE. Instead of SIS, we treat the bulge as a point 
mass, its value M e g is so selected for each q r that the predicted lensing probability can match 
the results of JVAS/CLASS represented by histogram. In the left panel, the value of the 
concentration parameter c\ is taken to be its median for any given halo mass and redshift. 
In the right panel, the scatter in c\ is considered by averaging the probability with the well 
known log-normal distribution. In both panels, og = 0.95, and the null result for lenses with 
6" < A# < 15" of JVAS/CLASS is shown with the thick dashed horizontal line indicating 
the upper limit. 
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log- normal distribution (right panel), the probabilities are increased considerably at larger 
image separations and only slightly increased at smaller separations for all cases, and the 
"scatter" among the four cases is reduced. So the scatter in c\ should be considered when 
one uses the NFW profile to constrain some related parameters. 

It is also shown in Fig. 1 that, for low flux ratios (q r < 10, as for the JVAS/CLASS 
survey), the NFW plus a single SMBH model produces far too few small separation lenses 
(the dash-dot-dot-dot line in the figure). This is confirmed by the fact that, among the 22 
confirmed lenses in JVAS/CLAS, none has a fainter image very close to the center of the 
lens (Browne et al. 2002). In other words, up to date strong gravitational lensing effect of a 
single SMBH has not been found. 

A larger flux ratio requires a smaller fraction of the bulge mass as a point mass. It's 
interesting to note that, when q r ps 10 4 and M c r ps 30M # , and when no scatter in C\ is 
included (the dotted line in the left panel of Fig. 1), the predicted lensing probability fit quite 
well the JVAS/CLASS results. As pointed out earlier, this conclusion can be equivalently 
obtained with the model of two populations of halos, the combination of SIS and NFW, when 
no constraints on flux density ratio q T are taken into account (i.e., q T ~ oo. Note that, the 
flux density ratios with values as high as q r ~ 10 4 and infinity have approximately the same 
effect on the predicted probabilities). However, since the flux density ratio for the well-defined 
JVAS/CLASS sample is limited to q r < 10, the above mentioned good fit cannot be regarded 
as reflecting the real nature of CDM halos. In fact, to match the JVAS/CLASS results, for 
q r < 10, about 20% of the bulge mass is required as a point mass (i.e., M e g ps 200M., the 
solid lines in Fig. 1), while for q r < 10 4 , only 3% of bulge mass is needed. This difference 
is very important when one attempts to use the observational results of JVAS/CLASS to 
constrain the density profile of galactic bulge and/or dark matter halos. It's well known 
that, the flux density ratio will reduce the lensing cross section considerably. All the models 
without considering this would have overestimated the lensing probabilities. In our NFW 
+ Point mass model, we can adjust the value of the point mass to the corresponding flux 
density ratio to match the observational results, and will never fail. In the NFW+SIS model, 
however, the contributions of SIS to lensing probabilities cannot be changed. Put in another 
way, after considering the flux density ratio limited by JVAS/CLASS, NFW+SIS model 
cannot reproduce the observational results, not even marginally, if some other parameters 
like M c are the same as used here. We have noted that, the solid lines in both panels of Fig. 
1 can only marginally fit the observations of JVAS/CLASS at larger image separations. But 
it is still acceptable, in the sense that, the curve is lower than the upper limit provided by 
the JVAS/CLASS survey for 6" < A9 < 15" (Phillips et al. 2001). We know that among 
the total of 22 confirmed 
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Fig. 2. — The same as the right panel of Fig. 1 except as for each panel here. From left to 
right, cr 8 is 0.7, 0.8, 0.9, 1.0 and 1.1, respectively. 
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lens systems of the JVAS/CLASS survey, 21 of them have image separations between 
0.3" and 3 ", and one of them (CLASS B2108+213) has an image separation of 4.6" (Browne 
et al. 2002). The null result of the JVAS/CLASS survey for 6" < A6 < 15" means that the 
NFW profile is suitable for all mass sizes of halos, when a certain fraction of the mass of 
each galactic bulge is considered. 

Although we believe that the reasonable match to the observations of JVAS/CLASS is 
q r < 10 and M e g/M, = 200 (the solid line in each panel of Fig. 1), it is helpful if we treat the 
four cases (four upper lines matching the histogram) as a whole to constrain the sensitive 
parameter ag. We plot in each panel of Fig. 2 the averaged lensing probability as a function 
of image separation A9. All the parameters are the same as those in the right panel of Fig. 
1 except for a s , for which five different values within the entire observational range (from 
0.7 to 1.1, as explicitly indicated) have been investigated to see their effect on the predicted 
lensing probabilities. Clearly, the larger values of a 8 will produce the higher probabilities. In 
the case of as = 0.7, too few lenses are produced at small image separations, so a low value 
of c"8 (< 0.7) is unlikely. The best fit is as ~ 1.0. This result is very close to that obtained 
most recently by Bahcall & Bode (2002), in which a 8 is determined from the abundance of 
massive clusters at redshifts z = 0.5 — 0.8. Our result is also in excellent agreement with 
that of Komatsu & Seljak (2002) (cr 8 = 1.04 ± 0.12 at 95%) suggested by the excess cosmic 
microwave background fluctuations detected on small scales by the CBI (Mason et al. 2002) 
and the BIMA (Dawson et al. 2002) experiments. 

In summary, strong gravitational lensing probability is quite sensitive to the cooling 
mass M c , concentration parameter Ci, flux density ratio q r and the spectrum normalization 
parameter as- Including the scatter in c\ in the calculations would increase the lensing 
probabilities considerably at larger image separations, while neglecting q r would overestimate 
the lensing probabilities. In our NFW + point mass model, 20% of the galactic bulge mass 
is required as a point mass to match the JVAS/CLASS results; the low values of a s (< 0.7) 
are ruled out, and the preferred value in our model is a 8 ~ 1.0. 

I thank Li-Xin Li for helpful discussions and the anonymous referee for useful sugges- 
tions. This work was supported by the National Natural Science Foundation of China under 
grant No. 10233040. 
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